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ABSTRACT 

Galaxy-galaxy interactions rearrange the baryons in galaxies and trigger substantial star formation; the aggre- 
gate effects of these interactions on the evolutionary histories of galaxies in the Universe are poorly understood. 
We combine B and R -band photometry and optical spectroscopy to estimate the strengths and timescales of bursts 
of triggered star formation in the centers of 190 galaxies in pairs and compact groups. Based on an analysis of 
the measured colors and EW(Ha), we characterize the pre-existing and triggered populations separately. The 
best-fitting burst scenarios assume stronger reddening corrections for line emission than for the continuum and 
continuous star formation lasting for > a hundred Myr. The most realistic scenarios require an initial mass function 
that is deficient in the highest-mass stars. The color of the pre-existing stellar population is the most significant 
source of uncertainty. 

Triggered star formation contributes substantially (probably > 50%) to the /?-band flux in the central regions of 
several galaxies; tidal tails do not necessarily accompany this star formation. Many of the galaxies in our sample 
have bluer centers than outskirts, suggesting that pre- or non-merger interactions may lead to evolution along the 
Hubble sequence. These objects would appear blue and compact at higher redshifts; the older, redder outskirts of 
the disks would be difficult to detect. Our data indicate that galaxies with larger separations on the sky contain 
weaker, and probably older, bursts of star formation on average. However, confirmation of these trends requires 
further constraints on the colors of the older stellar populations and on the reddening for individual galaxies. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: kinematics and 
dynamics — galaxies: interactions — galaxies: structure 



1. INTRODUCTION 

Larson & Tinsley (1978) showed that galaxies in pairs exhibit 
a broader scatter in the U - BIB - V plane than "field" galax- 
ies, demonstrating that interactions trigger bursts of star forma- 
tion. Since then, numerous studies, both of statistical samples 
and of individual interacting systems, have confirmed their re- 
sults (e.g., Kennicutt & Keel 1984; Madore 1986; Kennicutt et 
al. 1987; Jones & Stein 1989; Sekiguchi & Wolstencroft 1992; 
Keel 1993; Liu & Kennicutt 1995a,b; Keel 1996; Donzelli & 
Pastoriza 1997; Barton, Geller, & Kenyon 2000a). These stud- 
ies have conclusively established the link between interactions 
and star formation. However, many questions remain about the 
true role of interactions in the evolutionary history of galaxies 
in the Universe. 

Because interaction timescales are much shorter than a Hub- 
ble time, galaxies may participate in many interactions and/or 
mergers during their lifetimes. The strengths and durations of 
typical triggered bursts of star formation are therefore impor- 
tant cosmological parameters. An understanding of interactions 
is a necessary basis for including interactions and mergers in 
models of galaxy formation (e.g., Somerville & Primack 1999; 
Diaferio et al. 1999; Kauffmann et al. 1999a,b). Studies of lo- 
cal interactions also aid interpretation of the unusual morpholo- 
gies and structural parameters of intermediate redshift galax- 
ies (e.g., Koo et al. 1994; Abraham et al. 1996; Barton & van 
Zee 2001). Finally, a measure of the total effects of interac- 
tion may yield constraints on the causes of evolution along the 
Hubble sequence; large amounts of gas funneled into the cen- 
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ters of the galaxies could result in interaction-induced forms of 
the processes in secular evolution (Pfenniger & Norman 1990) 
that cause the formation of "exponential bulges" (Andredakis 
& Sanders 1994; Carollo 1999). 

Barton, Geller, & Kenyon (2000a, BGK hereafter) study 502 
galaxies in pairs or compact groups selected from the Cf A2 red- 
shift survey. They find a significant correlation between pah- 
separation on the sky, AD, and Ha equivalent width, EW(Ha). 
The correlation also extends to the AV/EW(Ha) plane, where 
AV is the pair separation in velocity. BGK argue that the AD - 
EW(Ha) correlation results from the aging of a continuing burst 
of star formation. If their interpretation is correct, the correla- 
tion provides a method of measuring the durations and initial 
mass functions (IMFs) of the bursts by comparing dynamical 
timescales with star-formation timescales. Here we investigate 
the origin of this correlation by using B and R photometry to 
constrain the old stellar population and the new burst of star 
formation separately. 

To measure the amount and current age of the new burst of 
star formation, we explore a method of using the measured 
colors and EW(Ha) to characterize recent star formation su- 
perposed on a pre-existing stellar population. We relate these 
quantities to the orbital parameters of the pairs, thus exploring 
the origin of the BGK correlation. In Sec. 2 we describe the 
pair sample and the data. Sec. 3 contains a description of the 
"two-population" model we apply to characterize the burst of 
star formation independent of the underlying older stellar pop- 
ulation. We include a brief discussion of the results of applying 
the model to the data. We discuss the origin of the BGK AD - 
EW(Ha) correlation in Sec. 4. In Sec. 5 we explore the depen- 
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dence of burst strength and age on the R-band luminosities and 
rotation speeds of the galaxies. We describe the galaxies with 
strong bursts of central star formation in Sec. 6 and conclude in 
Sec. 7. 

2. THE PAIR SAMPLE 

We base our study on the sample of all 786 galaxies in pairs 
and compact groups drawn from the original CfA2 redshift sur- 
vey (m Zw < 15.5) with AD < 50 h" 1 kpc, AV < 1000 km/s and 
v < 2300 km/s, where AV is the pair (or compact group neigh- 
bor) velocity separation, and v = cz is the apparent recession 
velocity. BGK estimate that this "full" sample of 786 galaxies 
from the original CfA2 survey is 70% complete with respect 
to all known galaxies in pairs over the region of interest in the 
updated CfA2 redshift survey (Falco et al. 1999). 41 of the 502 
galaxies with (new) spectra fail to satisfy the original selection 
criteria. We include these galaxies in the analysis because most 
come close to satisfying the criteria; all pairs in the sample sat- 
isfy AD < 77 h" 1 kpc and AV < 1035 km/s. 

We divide the full sample into three overlapping sub-samples 
for which we have different kinds of data: the dynamical sam- 
ple (~140 galaxies, mostly spirals and SOs, with Ha rota- 
tion curves), the spectroscopic sample (502 galaxies with one- 
dimensional spectra; the same as the BGK sample) and the pho- 
tometric sample (196 galaxies with calibrated photometric im- 
ages in B and R). 

In this paper, we concentrate on the photometric properties 
of the galaxies in the overlap of the spectroscopic and photo- 
metric samples (190 galaxies comprising all or part of 78 pairs 
and 13 compact groups; SP sample hereafter). We have Ha ro- 
tation curves for 103 of these galaxies; 89 provide useful mea- 
surements of V c , the velocity width. BGK separate the sample 
into galaxies residing in low- and high-density environments 
and concentrate on the low-density subset to avoid contamina- 
tion by the morphology-density relation and gas stripping in 
clusters. Here our sub-samples are much smaller and we are 
unable to divide the data in this way. However, we use direct 
measurements of V c to understand the galaxy mass-dependent 
effects. 

In spite of the m Zw < 15.5 cutoff, the pair sample contains 
numerous galaxies with B magnitudes > 15.5. These systems 
were often recorded in the CGCG (Zwicky & Kowal 1968), and 
therefore the redshift catalog, with the Zwicky magnitude of the 
whole system instead of separate magnitudes for each of the 
galaxies. We find 33 (of 190) galaxies (17%) with m B 26 > 15.5 
in the SP sample, where m B 26 is the fi s = 26 isophotal magni- 
tude. The sample is incomplete with respect to these galaxies. 

We selected the systems in the dynamical sub-sample based 
on the suitability of at least one galaxy for a rotation curve mea- 
surement; therefore this sub-sample heavily favors Ha-emitting 
galaxies and spirals less than ^50° from edge on. For the selec- 
tion of the photometric and spectroscopic samples, we focus on 
(1) galaxies in the dynamical sample, and (2) galaxies known 
to have emission lines. Therefore, the SP sample heavily favors 
emission-line galaxies. We measure the extent of the bias by 
comparing the spectroscopic properties of the SP sample with 
those of the CfA2South pair sample, which is 100% complete 
with respect to the set of all pairs selected from the original 
CfA2 survey. Although only 25% of the complete CfA2South 
sample has detectable EW(Ha), 75% of the (overlapping) SP 
sample has detectable EW(Ha). 

The SP sample does not favor line-emitting galaxies with 
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FIG. 1. — The relationship between AD and EW(Ha) for the 190-galaxy 
SP sample (for comparison with BGK Fig. 2a). The trend is significant, with 
P SR = 0.058 in the full SP sample and P SR = 0.0017 in the sample restricted 
to galaxies with measured Balmer decrements (filled circles). We exclude 2 
galaxies with AD > 55 h _1 kpc from the plot. 



particularly strong emission. The distribution of EW(Ha) in 
the line-emitting fraction of the SP sample is unbiased with re- 
spect to the EW(Ha) distribution for the complete CfA2South 
sample; a K-S test is unable to detect differences in the EW(Ha) 
distribution of galaxies with EW(Ha) > 10 A in the (overlap- 
ping) SP sample (115 galaxies) and the CfA2South spectro- 
scopic samples (98 galaxies, with 44 in the SP sample), with 
Pks = 0.55 including all galaxies, or Pks = 0.18 excluding the 
overlapping galaxies. We plot AD vs. EW(Ha) (uncorrected 
for reddening) for the SP sample in Fig. 1. 

Differences are evident in the AD and redshift distributions 
of the SP and complete CfA2South samples. For example, the 
SP sample contains an excess of pairs with AD ~ 10 h" 1 kpc 
and an excess of pairs at lower redshift. These effects result 
from the large-scale structure differences between CfA2North 
and CfA2South (see Falco et al. 1999). We use Spearman rank 
tests to evaluate the significance of the correlations we describe 
below. Therefore, the differences in the redshift and AD dis- 
tributions should only have a significant bearing on our conclu- 
sions if we have different selection effects in different ranges of 
AD. We are not aware of any such effects. 

2.1. Observations and Data Reduction 

BGK and Barton et al. (2001) describe the basic observations 
and data reduction; we summarize these observations here and 
we list the results in Table 1 . 

We observed the galaxies at the FLWO 48" telescope on Mt. 
Hopkins through either Johnson or Harris B and Cousins or 
Harris R filters, for total exposure times of 15 minutes (usu- 
ally spread over 3 images) and 5 minutes (usually exposed all 
at once or spread over 2 images), respectively. After calibrating 
the images to the Landolt (1992) standard-star system, we sub- 
tract the sky background and measure magnitudes and surface 
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Table 1 
The SP Sample 



Coordinates 




cz 


AD 


Ay 






EW(Ha) 




v 2.2 




(1950) 


N S y S 


(km/s) 


(kpc/h) 


(km/s) 


Mr 


(B-R) s 


(A) 


Ha 

H<3 


(km/s) 


Note Label 


00 03 47.2+17 09 05 


2 


5590 


7.6 


181 


-19.94* 


0.75 


130.1 


4.07 




marg 


00 03 47.5 +17 09 32 


2 


5770 


7.6 


181 


-20.14* 


0.92 


47.1 


4.82 


221.5 


norm 


00 1 1 07.7 +47 52 24 


2 


5445 


31.9 


22 


-21.40* 


1.39 


0.0 








00 11 19.5+47 52 47 


2 


5422 


31.9 


22 


-21.17* 


1.51 


5.8 




372.6 


norm 


00 17 52.0-00 32 55 


2 


5222 


9.8 


423 


-20.84* 


1.38 


9.0 


2.88 


247.9 


marg 



"Image was calibrated from a shorter exposure. 
^Galaxies overlapped. 

c Bright star in background, or smaller but significant star(s) on galaxy. 

Note. — Galaxies in the SP sample: (1) 1950 coordinates, (2) number of galaxies in system, (3) redshift, (4) separation of pair on the sky (or distance to closest 
galaxy in n-tuple), (5) velocity separation of pair (or velocity difference to closest galaxy in velocity within the n-tuple), (6) R absolute magnitude, extrapolated to 
the total magnitude and measured assuming Ho = 70 km s -1 Mpc~' , (7) color of the portion of the galaxy on the slit during the spectroscopic observation, corrected 
for Galactic extinction and nebular lines, (8) Ha equivalent width, (9) Balmer decrement (corrected for Balmer absorption), (10) corrected velocity width, and (11) a 
description of rotation curve, where norm = "normal", marg = "marginal" (intermediate between normal and distorted), and dist = "distorted" (see text of Barton et al. 
2001), and (12) label for Figs. 23 and 24. [The complete version of this table will appear in the electronic edition of the Journal and is available from the first author.] 



brightness profiles with ellipse fitting on the B-band images; we 
overlay B-band isophotes onto the B-band images. Barton et al. 
(2001) discuss the Galactic reddening correction and several is- 
sues relating to galaxy distortion and overlap. 

BGK contains a complete description of the spectroscopy 
and data reduction. We observed the galaxies using the FAST 
spectrograph at the 1.5m Tillinghast telescope on Mt. Hopkins. 
Typical exposures lasted ~10 - 20 minutes. We observed the 
galaxies through a 3" slit using a grating with 300 lines/mm to 
disperse the light between 4000 and 7000 A. The galaxies were 
often observed through clouds; thus, the spectra are roughly 
calibrated only in a relative sense, using observations of one or 
more spectrophotometric standards each night. The photometry 
of the galaxies provides a more reliable indicator of the over- 
all continuum shape of the spectral energy distribution. Fig. 2 
shows examples of the spectra, illustrating the observed range 
of star-forming properties. 

After fiatfielding the data, we extract apertures which in- 
clude the brightest continuum light incident on the slit from the 
galaxy. This procedure results in a range of aperture lengths, 
from 0.25 - 13.7 h" 1 kpc, with a median of only 2.4 h" 1 kpc. 
Thus, strictly speaking, the spectra are not "nuclear" but "cen- 
tral" spectra encompassing the central light distribution. We 
compute equivalent widths that are the ratio of the flux in the 
line to the measured continuum. We base errors on photon 
statistics; repeat measurements indicate that these errors are 
within a factor of ^2 of the true measurement errors. We cor- 
rect for Balmer absorption by computing the maximum absorp- 
tion equivalent width in HS or H7 and adding that equivalent 
width to both Ha and H(3. 

For the SP sample, we use the cleaned 2 48-inch images (see 
Barton et al. 2001) to measure the B and R isophotal mag- 
nitudes, which we extrapolate to total magnitudes, and B-R 
color, B s , R s , and (B—R) s , respectively, of the light that accumu- 
lated through the slit during the spectroscopic observations. We 
compute these quantities by adding the flux within a rectangular 
aperture measuring 3" x L np and aligned with the spectroscopic 
aperture (to within < 1 .5°), for both the B and R clean images. 
3" is the width of the slit and L ap is the (variable) length of the 
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FIG. 2. — Sample spectra illustrating the range of observed star-forming 
properties of the interacting galaxies. To show the quality of the data, the 
spectra have not been smoothed. 



extracted aperture. We assume the slit center and the center of 
the galaxy (the center of the smallest fit ellipse) are coincident, 
a close approximation to the actual position of the slit during 
the spectroscopic observation. We measure conservative errors 
in (B-R) s due to misplacement of the slit center by displacing 
the slit 2 pixels (~ 1.2") in the ±x and ±y directions, for a 
subset of 30 galaxies. Changes in the color range from ^-0. 1 to 
^+0.2 magnitudes, with a = 0.055. When the spectrum consists 
of two or more exposures, usually with different values of L ap , 
we scale the fluxes in the different apertures by the exposure 
times of the spectra; this procedure is only approximately cor- 
rect in cases where the atmospheric absorption differed during 
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the observations. 

Finally, we observed galaxies in the dynamical sample with 
the Blue Channel Spectrograph at the Multiple Mirror Tele- 
scope on Mt. Hopkins between November 1996 and February 
1998. For the majority of the galaxies, we used a 1" slit and 
a grating with with 1200 lines/mm centered at ^6500 A, and 
including roughly the wavelength range 5800 - 7200 A. 

Barton et al. (2000b) and Barton et al. (2001) describe the 
data reduction in detail. We use a cross-correlation technique 
that makes simultaneous use of the major emission lines: [Nil] 
(A6548 and A6583), [SII] (A6716 and A6730), and Ha; this 
technique is straightforward to implement and yields well- 
defined errors. We use the method of Courteau (1997), de- 
veloped for Tully-Fisher (Tully & Fisher 1977, TF hereafter) 
applications, to measure the velocity width, V2.2, for galaxies 
in the dynamical sample. We fit each rotation curve with an 
empirical fitting function (see Courteau 1997, Eqn. 2). Barton 
et al. (2001) describe the corrections we apply for inclination, 
slit misalignment, and redshift, to arrive at N c 2 2 , the corrected 
velocity width. 

3. PARAMETERIZING THE CENTRAL BURSTS OF STAR 
FORMATION 

For lack of adequate constraints, full modeling of an unre- 
solved stellar population requires assumptions about the star 
formation rate as a function of time. Many studies parameterize 
the star formation rate with a single function, such as an expo- 
nential, and fit for the timescale, r (e.g., Larson & Tinsley 1978; 
Abraham et al. 1996; Bell & de Jong 2000); others find the best- 
fit spectral energy distribution from a set of empirical templates 
(e.g., Thompson, Weymann, & Storrie-Lombardi 2001). Be- 
cause our goal is to parameterize a recent (< 1 Gyr) burst of star 
formation superposed on an otherwise typical galaxy, we must 
parameterize the "old" (pre-existing) and "new" (interaction- 
related) populations separately. 

Several existing studies address this "added-burst" problem 
(e.g., Larson & Tinsley 1978; Kennicutt et al. 1987; Kennicutt, 
Tamblyn, & Congdon 1994). For all of these studies, there is 
a degeneracy associated with separating the new burst from the 
pre-existing galaxy: is the galaxy blue because its star forma- 
tion rate has been constant or increasing with time throughout 
its history, or is it blue because of a strong triggered burst? Or, 
cast in terms of the pair population as a whole, does the range 
of colors observed in pairs result from the range of colors in the 
progenitor galaxies or from triggered bursts of star formation? 
One approach involves beginning with the distribution of col- 
ors vs. EW(Ha) observed in the field and exploring the effects 
of short bursts of star formation (e.g., Kennicutt et al. 1987). 
However, every point in the color/EW(Ha) plane is reachable 
from a range of starting points. In an approach similar to the 
models of Larson & Tinsley (1978), we essentially fix assump- 
tions about the pre-existing stellar populations and characterize 
bursts of star formation in terms of burst age and "strength", 
testing the effects of modifying the initial assumptions. We 
avoid assumptions about mass-to-light ratios by defining burst 
strength in terms of the fraction of the /?-band flux originating 
from the new burst. In Sees. 3.2-4 we discuss the full range of 
possibilities; we include a discussion of techniques for resolv- 
ing the degeneracy between young pre-existing stellar popula- 

2 The "cleaned" images have been stripped of overlapping stars and, in some 
cases, overlapping companion galaxies; we usually replace the stripped regions 
with an elliptically symmetric model of the galaxy. 



tions and triggered bursts. We also discuss the effects of other 
properties, like metallicity, on our results. We demonstrate the 
effectiveness of using only B-R and EW(Ha) in Sec. 3.3. 

3.1. The Two-Population Model 

For simplicity, we begin by assuming that two distinct stel- 
lar populations contribute to the flux incident on the slit aper- 
ture during the spectroscopic observation: (1) a new popula- 
tion, solely responsible for any Ha emission, which we denote 
by the subscript n, with B and R fluxes, fB n (t) and fR„(t), respec- 
tively, that turn on at t = and evolve as the burst ages, and (2) a 
pre-existing, older population, denoted by the subscript 0, with 
fluxes fso and fRo, which do not evolve significantly on the burst 
timescales (< 1 Gyr). We explore the effects of Ha emission 
from the pre-existing stellar population in Sec. 3.2. With these 
assumptions, the Ha equivalent width from spectroscopic mea- 
surements [EW(Ha) mea s(t)], slit B and R fluxes from calibrated 
photometry, and starburst models for simple burst histories, we 
estimate both the strengths and ages of the starbursts. 

The old population increases the continuum around Ha. 
Therefore, the fractional strength of the (new) burst in R, Sr(0 = 
ho+hV) ' g° verns trie relationship between EW(Ha) meas (t) and 
the Ha equivalent width of the burst alone, EW(Ha)b urs t(t)- The 
measured equivalent width is 

EW(Ha) meas (t) = s R (t) x k(z,SED ,SED b ) x EW(Ha) burst (t), 

(1) 

where k(z, SEDo, SEDb) is a factor of order unity (0.9 < 
k(z,SEDo,SEDh) < 1.4) that corrects for the difference be- 
tween the burst strength in R and the burst strength at the con- 
tinuum near Ha. k(z, SEDo, SEDb) depends on the spectral en- 
ergy distributions of the old stellar population (SEDo) and the 
burst (SEDb). We estimate k(z, SEDo, SEDb) at the appropri- 
ate redshift for a given old stellar population and burst by using 
Bruzual & Chariot (1996) Salpeter r-model spectral energy dis- 
tributions for the old population, matched in B-R to the nearest 
t model at 14 Gyr from formation (with r equal to an integral 
number of Gyr in most cases). We interpolate the Leitherer et 
al. (1999) burst model spectral energy distributions for the burst 
SED computed at a limited set of burst ages. We use the R filter 
function of Bessell (1990) and the B filter function of Buser & 
Kurucz (1978). In effect, we parameterize the old stellar pop- 
ulation completely via (B - R)o which is an input parameter, 
and k(z, SEDq, SEDb), which we compute from the Bruzual & 
Chariot (1996) models based on the input (B-R) Q . 

The metallicity of the old stellar population has little effect 
on k(z, SEDo, SEDb). For a fixed (B-R)o = 1.5, we estimate 
the effects of metallicity by computing k(z,SED ,SED b ) over 
the full range of strengths from to 1 and the full range of 
ages, IMFs, metallicities, and star foramtion histories of the 
L99 models. The maximum difference in k(z, SEDo, SEDb) 
introduced by assuming a non-solar metallicity (of either 0.4 
Z or 2.5 Z ) is < 1%. Hereafter, we use Z= Z© to compute 
k(z, SED , SED b ). 

The strength of the new burst also depends on the colors of 
the old and new stellar populations. The composite (measured) 
R-to-B flux ratio is 

Cs(t) = 1 1 -, M 'Co-CnWT ' (2) 

l + s R (t) Cn(t) 

where Co and C n (t) are the R-to-B flux ratios of the light in- 
cident on the slit from the old and new populations respec- 
tively. The R-to-B flux ratios are related to the colors: Co = 
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Table 2 

Assumptions for Burst Model Contours 



Parameter 


Assumption 
for Fig. 3 


Other assumptions 
from this paper or BGK. 


Pre-existing population 
color, (B-R) 


1.5 


0.9-2.0 
1.5 + 0.75 x [log 2 -2.5] 
1.5 + 1.50 x [log Y° 2 2 -2.5] 
Peletier & Balcells (1996) bulges 
NFGS distributions 


Burst IMF 


"steep" (a = 3.3) 


"Salpeter" (a = 2.35) 
"cutoff Salpeter" (a = 2.35, Mass < 30 M Q ) 


Burst SFR(t) 


Constant for 1 Gyr 


Instantaneous 


Burst metallicity 


Solar 


O.4Z ,2Z Q 



lm = 10 (^)o/2.5 ) Cs(t) = (|2±M2) = lO^V 2 - 5 , and C n (t) = 

/Rn(t) _ 1 Q(S-K) n /2.5 
An(t) 

The "strength" parameter in Eqns. 1 and 2, s R (t), is time- 
dependent. It represents the current fraction of the central R- 
band flux originating from the new burst; this fraction increases 
as the burst builds luminosity and decreases as the burst fades. 
We also consider a time-independent parameter, sioo, the flux 
at an age of 100 Myr [s 100 = s R (100Myr)]. Then at any other 
time, 

_Sioo_ jq0.4[Mr(100 Myr)-M R (t)] 

„ /|\ _ 1~ S 1Q0 /o\ 

SRW - j + sioo jq0.4[M r (100 Myr)-M R (t)] ' ^ ' 

1— SlOO 

where M R (t) is the /?-band absolute magnitude of the burst as a 
function of time. 

3.1.1. Strengths and Ages from the Leitherer et al. Models 

Starburst models predict the burst color, (B-R) n , and the 
burst Ha equivalent width, EW(Ha)burst(t), as a function of 
time. We use the Leitherer et al. (1999, L99 hereafter) starburst 
models for burst parameters. L99 model a starburst by "form- 
ing" stars according to an IMF with |g = CM"", where M is 
the stellar mass, a is the slope of the IMF, and C is a constant. 
With models from the literature or their own models for stel- 
lar evolution, stellar atmospheres, and nebular emission, they 
compute EW(Ha) as a function of time for both the "instanta- 
neous" starburst case, with a timescale t < 10 6 years, and the 
continuous case, with star formation over 10 9 years. 

The L99 models have a range of metallicities (0.05 Z©, 0.2 
Z©, 0.4 Z©, Z©, 2Z Q ), and a range of IMFs, including an 
approximate Salpeter (1955) IMF, with a = 2.35, from 1 - 
100 M©, a Salpeter IMF with a 30 M© high-mass cutoff, and 
a steeper function with a = 3.3 from 1-100 M© patterned af- 
ter the high-mass segment of the Miller & Scalo (1979) IMF. 
Hereafter, we refer to these IMFs as Salpeter, cutoff Salpeter, 
and steep, respectively. 

For an assumed pre-existing stellar population, the starburst 
models and Eqns. 1-3 provide a method of estimating burst 
strengths, s R (t) or sioo, and ages, t, from the measured quanti- 
ties [(B-R) s and EW(Ha) meas (t)]. Fig. 3 shows contours of con- 
stant current burst strength, s R (t), and age, t, in the (measured) 
EW(Ha)/(B-fi) s plane for a continuous burst of star formation 
with the "steep" IMF and (B-R)o =1.5. Hereafter, we assume 
cz = 5000 km s" 1 to compute k(z, SEDo, SEDb), except for spe- 



cific galaxies, for which we use the measured redshift. Table 2 
summarizes the other assumptions we make to construct these 
contours. For comparison, Fig. 4 shows the contours of con- 
stant (current) burst strength assuming old population colors of 
(B-/?)o = 1.0and(B-/?) = 2.0. 
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FIG. 3. — Contours of constant current burst strength and age in the plane of 
the observed parameters, for a continuous burst of star formation (see Table 2). 
The dotted lines are lines of constant burst age, t, in log(years); solid lines are 
contours of constant current strength, SR(t), defined as the fraction of R-band 
flux from the new burst in the spectroscopic aperture. 

Lines of constant tb urs t are diagonal in Figs. 3 and 4, while 
the contours of constant s R (t) are nearly vertical. Thus, for con- 
tinuous bursts of star formation with steep IMFs superposed on 
fixed existing stellar populations, the burst age, tburst, primarily 
affects the measured EW(Ha). In contrast, the strength of the 
burst in R, s R (t), affects both the measured EW(Ha) and the 
measured color, (B-R) s . 

Contours of constant sioo in Fig. 5 are the evolutionary tracks 
of constant bursts of star formation within this model, sioo is 
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FIG. 4. — Same as Fig. 3, but with (B—R)q = 1.0 (left) and (B—R)q = 2.0 (right). Lines of constant current strength, srO:), remain roughly vertical, but the lines 
corresponding to weaker bursts are displaced depending on (B-R)q. The blue edge of the contours is independent of (B-R)q; the IMF and the starburst models 
determine the location of a population that is 100% triggered burst. 



the burst strength in R that would be measured at 100 Myr, as- 
suming that the burst continues (or has continued) at a constant 
rate. Sioo depends heavily the assumption of a constant star for- 
mation rate. Mass-fraction measures (e.g., Larson & Tinsley 
1978; Kennicutt et al. 1987) require similar assumptions; they 
incorporate assumptions about the relative mass-to-light ratios 
of the burst and the underlying old stellar population. 

The contours in Fig. 5 fold over and are no longer verti- 
cal. A moderately strong, young burst continuing at its current 
rate dominates the flux of the galaxy at 100 Myr. For exam- 
ple, within this model, if a continuous burst superposed on an 
older population [(B-R) Q = 1.5] has EW(Ha) meas (t) = 120 A 
and (B-R) s = 0.6, its current age (Fig. 3) is ^16 Myr and ^50% 
of the /?-band flux currently originates from this burst. If this 
burst continues at its present rate, at 100 Myr it will comprise 
^80% of the /?-band flux in the aperture (Fig. 5). If the burst 
actually ceases at 316 Myr, it will comprise much less than the 
original 50% of the R-band flux in the aperture at 100 Myr. 

3.2. Measuring Strengths and Ages from the Data 

The discussion in Sec. 3.1.1 omits several corrections nec- 
essary to apply the models to the data. Table 3 summarizes 
the corrections we discuss in this section. They include (1) 
the contribution from line emission to (B-R)^ which we esti- 
mate using the measured equivalent widths of the nebular lines 
along with the continuum from the model spectral energy dis- 
tributions. We apply the nebular emission (and absorption) cor- 
rections to the color data in the slit regions only (not to the 
total magnitudes), and include the corrections in Figs. 3-5; 
the corrections range from -0.079 to 0.003 with an average 
of -0.017, (2) extinction corrections to the measured colors 
and/or EW(Ha) meas (t), and (3) ongoing star formation in the 
pre-existing stellar population. Finally, we discuss our choices 
for the color of the pre-existing stellar population, (B- R)o- 



Table 3 
Corrections to the Data 

Correction(s) explored 
Correction in this paper 

nebular line contamination in color corrected data using measured 

equivalent widths, model SEDs 
Internal reddening "traditional'" 1 

none 

"Calzetti" b 

ongoing star formation in no correction applied 

pre-existing stellar population 



"See Sec. 3.1. 

^Calzetti, Kinney, & Storchi-Bergmann (1994); see Sec. 3.1. 



3.2.1. Corrections for Extinction 

Corrections for internal extinction substantially modify the 
inferred ages and strengths of the bursts of star formation. Be- 
cause it is independent of the star formation history of the 
galaxy, the Balmer decrement 3 is generally the best extinction 
measure available in our dataset. However, several circum- 
stances affect our ability to measure the Balmer decrement and 
to derive the appropriate extinction correction for the galaxy. 
The presence of an intermediate-age (^Gyr) population neces- 
sitates a correction for Balmer absorption under the emission 
lines (see Sec. 2. 1 and BGK for a description of our correction). 
In addition, Calzetti, Kinney, & Storchi-Bergmann (1994) show 
that the extinction of the continuum is only half of that appro- 
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priate for line-emitting regions. 

The difficulty of measuring the Balmer decrement in galax- 
ies with insufficient star formation introduces selection effects 
into the extinction-corrected sample. The top panel of Fig. 6 
shows the fractional error in the Balmer decrement as a func- 
tion of EW(Ha) for the 78 galaxies with accurate measure- 
ments (better than 25%); 112/190 galaxies yield larger errors 
in Ha/H/3 or exhibit no significant H/3 emission. The bot- 
tom panel shows a histogram of EW(Ha) meas (t) for the galaxies 
with no Ha/H/3 measurements; incompleteness in the sample 
of galaxies with accurate Balmer decrement measurements be- 
gins at EW(Ha) mea s(t) ~ 20 A. In the analysis below, we pri- 
marily use the 78 galaxies with accurate measures; we replace 
values of Ha/H/3 < 2.85 with 2.85 (or 2.88 as necessary). 

We explore the possible effects of dust on our data by con- 
sidering three different dust corrections: no correction, a "tra- 
ditional" correction with no differences in extinction between 
the continuum and the line-emitting regions, and the extinc- 
tion correction of Calzetti, Kinney, & Storchi-Bergmann (1994, 
"Calzetti" correction, hereafter). 

The "traditional" correction involves no adjustments to 
EW(Ha) meas(t)- We correct (B — R) s for reddening assuming 

E(B-R) = 1.78E(fi-V)= 1.536E(/3-a) 

"(Ha/H/3)" 



1.536 x 2.5 log 



2.85 



(Zombeck 1990; Miller & Mathews 1972; Veilleux et al. 1995). 
Here, E(B - R), E(B-V), and E(/3-a) are the B-to-R, B-to-V, 
and H/3-to-Ha reddening, respectively, and Ha/H/3 is the flux 
ratio of Ha to H/3, corrected for Balmer absorption. The error in 

3 The Balmer decrement is the flux ratio Ha/H/3, where both fluxes have 
been corrected for stellar absorption with the prescription described in Barton, 
Geller, & Kenyon (2000a). 



the measured Ha/H/3 includes a small addition for uncertainty 
in the Balmer absorption correction. 

The "Calzetti" correction is appropriate if dust is clumped 
around regions of recent star formation (e.g., Chariot & Fall 
2000); the prescription requires separate corrections for the 
continuum and line emission (Calzetti, Kinney, & Storchi- 
Bergmann 1994): 

"Ha/H/3" 



and 



E(B-R) = 0.8170 In 



EW(Ha) c = EW(Ha) meas (t) x 



2.88 



Ha/H/3 
2.88 



0.9005 



3.2.2. Star Formation in the Pre-existing Stellar Population 

In this section, we discuss our choice not to correct for on- 
going star formation in the pre-existing stellar population. The 
appropriate corrections depend on the intended interpretation of 
the results. The fundamental question is whether the new burst 
of star formation is physically independent of the ongoing star 
formation in the progenitor galaxy. If the pre-existing stellar 
population and the recently-formed burst population are truly 
independent of one another, ongoing star formation in the pre- 
existing ("old") population necessitates corrections to the data 
and models. Non-zero corrections to Ha result in weaker burst 
strengths and older burst ages. We argue, however, that once an 
interaction disrupts a galaxy, star formation is not independent 
in the pre-existing and burst populations. The process that trig- 
gers the new, added burst of star formation likely disrupts the 
"quiescent" star formation in the galaxy. The triggered burst 
may consume the gas that was forming stars in the absence of 
an interaction, invalidating the assumption of physically inde- 
pendent pre-existing and burst populations. 
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FIG. 5. — Same as Fig. 3, except that solid lines are contours of constant 
strength for the burst at 100 Myr. Thus, constant sioo lines track the evolution 
of a continuing burst of star formation. Unlike sr(i), the measurement of sioo 
depends heavily on the assumption the burst star formation rate is constant. 
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FIG. 6. — (Top) Fractional error in the Balmer decrement as a function of 
EW(Ha) meas (t). Our cutoff is an error of 25% of the measured value. (Bot- 
tom) Histogram of the galaxies with > 25% errors. The set of galaxies with 
measured Balmer decrements is incomplete for EW(Ha) mKls (t) < 20-30 A. 
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FIG. 7. — An example of the model spectrum and the galaxy spectral en- 
ergy distribution for the "clean" case of NGC 2719A, which has a Balmer 
decrement consistent with no reddening. The thick black line is the spectral 
energy distribution of the model, fit to the B - R color and EW(Ha) of the 
center of the galaxy only. The model is the sum of the relevant L99 model, a 
21-Gyr-old burst with a "steep" IMF, and the Bruzual & Chariot (1996) model 
with a Salpeter IMF and an exponential star formation timescale of r = 3 Gyr, 
observed 14 Gyr after formation. The other nebular emission lines are not 
modeled. We plot broad band fluxes as filled circles; R is normalized to the 
spectrum and B is plotted at its effective wavelength. The data and models 
have both been smoothed by 10 bins. 



Without physical distinctions between the pre-existing and 
triggered populations, the amount of triggered star formation 
is an ill-defined quantity. Two better-defined, distinct quanti- 
ties are: (1) the statistical increase in central star formation due 
to major interactions, and (2) the total star formation ongoing 
in the centers of interacting galaxies (including star formation 
that would have occurred in the absence of a major interaction). 
Measuring these quantities requires different approaches to the 
correction for ongoing star formation in the pre-existing stellar 
population. The first requires a statistical correction for quies- 
cent star formation; the second question requires no correction, 
but a judicious choice for (B-R)o that allows for the possibility 
of star formation just prior to the interaction. 

Regardless of the interpretation, a moderate amount of ongo- 
ing star formation in the pre-existing stellar population has few 
qualitative effects. The primary effect is to place galaxies with 
no interaction-triggered star formation onto the apparent burst 
contours (as in Fig. 3). These galaxies occupy a locus across 
the lower portion of the (B - /?) s /EW(Ha) meas (t) plane defined 
by "normal" galaxies. For these galaxies, the model analysis is 
not applicable. Thus, some of the galaxies that appear, in this 
model, to have weak bursts of star formation may actually be 
quiescent late-type galaxies; our analysis is more reliable for 
galaxies with larger EW(Ha) meas (t). 

In the sections that follow, we focus on the approach of char- 
acterizing the physical conditions in galaxies with triggered star 
formation. We apply no corrections for ongoing star-formation 
directly to the data, but explore the effects of different choices 
for (B -/?)(,. 
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3.3. The Scope of the Model 

To characterize the star-forming properties of our galaxies, 
our approach is to fit simple star-formation models to the ob- 
served B-R colors and EW(Ha) in the centers of the galaxies, 
using spectral synthesis models from the literature (Leitherer 
et al. 1999; Bruzual & Chariot 1996). These models predict 
the entire spectrum of the galaxy, with the exception of most 
of the nebular line strengths; the agreement between the pre- 
dicted and observed spectra serves as a consistency check of 
our models. Fig. 7 shows the "clean" case of NGC 2719A, 
with a measured Balmer decrement that yields no evidence for 
reddening in the spectrum. The thin solid line is the data; the 
thick solid line is sum of the model spectra from Bruzual & 
Chariot (1996) and L99, respectively; we also include the Ha 
line, with the appropriate measured EW(Ha). For the old stel- 
lar population {dashed line), we use a ? r = 3 model with a 
Salpeter IMF at 14 Gyr, corresponding to (B-R) = 1.46 [close 
to our fiducial (B-R) Q =1.5]. The measured (B-R) s = 0.56 
and EW(Ha) = 100 A yield a burst strength and age of 0.51 
and 21.4 Myr, respectively, for this (B — R)o and the "steep" 
IMF, with solar metallicity. Hence, the "new burst" spectrum 
(dotted line) is the L99 model with continuous star formation 
and "steep" IMF after 21 Myr. The thick black line is the sum 
of these models, combined according to the derived value of 
Sr(i). Filled circles represent measured B-R color, plotted at 
the effective wavelength of each band with both filter passbands 
from Bessell (1990); we normalize R and plot B. By design, the 
broad-band color agrees perfectly with the model. 

Apart from the nebular emission lines, which are not pre- 
dicted by the models, the spectral shape agrees well with the 
model, even though it is not used as an input. This agreement 
results from the agreement between the rough calibration of the 
spectrum and the more reliable broad band colors, which are 
generally within ^0.2 magnitudes of one another for the galax- 
ies with no reddening. The exception is the very blue end, be- 
low ^4000 A, where the calibration of the spectrum is highly 
uncertain. The good agreement between the predicted spec- 
trum and the observed spectrum justifies our use of only the 
broad band colors and EW(Ha) a posteriori: the spectra do 
not contain any more of the information predicted by the model 
than we are already using. In particular, the model spectra cor- 
responding to other (B - R)q values [and adjusted according to 
the new derived Sr(i) and tb urs t] appear to agree similarly well, 
within the limits of the expected accuracy of the calibration of 
the spectrum, indicating that the optical spectra contain no ad- 
ditional reliable information about the color of the pre-existing 
stellar population. 

Broadly speaking, the spectral energy distribution of a galaxy 
with a triggered burst of star formation is determined by a large 
number of factors, including the metallicities, star formation 
histories, and IMFs of both the old and new stars, and the red- 
dening appropriate for each stellar population. As indicated 
above, the current dataset does not provide separate informa- 
tion about all of these quantities. Our approach is to identify 
the most significant factors and focus on their effects. As we 
describe in Sec. 3.1, for the fiducial (B-R)q = 1.5, the metallic- 
ity of the old stellar population only affects k(z, SEDo, SEDb) by 
< 1%. In Sec. 3.4, we characterize the effects of the metallic- 
ity of the burst and show that they are dwarfed by uncertainties 
in the burst IMF and old stellar population color. Hence, we 
focus on the most important factors, the star formation histo- 
ries, burst IMF, and reddening corrections. Our approach, as 
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FIG. 8. — Solutions for burst strength, SR(t) (left), and age, tburst (right) for NGC 2719A as a function of (B-R)q, for the three different IMFs, assuming solar 
metallicity. NGC 2719A has EW(Ha) = 100 A and (B-i?) s = 0.56 with no evidence for extinction in the optical spectrum. 



described above, is to parameterize the burst star formation his- 
tory in terms of an age, tburst, and a current /?-band strength, 
Sr(i), and compute the results for a set of different assumptions 
about the burst IMF, the reddening, and the old stellar popula- 
tion star formation history and IMF [as characterized together 
by (B-R)o]. Our argument is statistical in nature; we attempt 
to account for the relationship between the ensemble properties 
of the sample and the predictions of the models. 

3.4. Results 

The measured s R (t) and tb urst depend relatively strongly on 
the input color of the old stellar population, (B-R)o. Fig. 8 
shows this dependence for NGC 27 19 A, the galaxy depicted 
in Fig. 7 with a strong burst of star formation and no apparent 
reddening. The realistic possibilities for (B-R)o do not include 
the entire range in the figure (from 0.5 - 2.5). For galaxies of 
types SO to Sbc, Peletier & Balcells (1996) find bulge colors 
in the range ofl.l5<fi— /?<1.8; centers of later-type spirals 
and irregulars are as blue as ~0.8 for "field" galaxies (Bell & 
de Jong 2000). 

The uncertainties in the assumptions about IMF and (B-R)o 
dwarf the uncertainties introduced by assuming solar metallic- 
ity for the old and new bursts. As discussed in Sees. 3.1 and 
3.3, uncertainties introduced by the metallicity of the old stel- 
lar population are negligible. Fig. 9 illustrates the typical ef- 
fects of our assumption of solar metallicity for the burst; the 
effects are small compared with the other uncertainties. To ver- 
ify that the observed range of metallicities in our sample is lim- 
ited, we use the "combined" metallicity estimator of Kewley & 
Dopita (2002); 99/190 of the SP sample spectra exhibit strong 
enough emission lines to make this measurement. For a so- 
lar abundance of \og(0/H)+ 12 = 8.69 (Allende Prieto, Lam- 



bert, & Asplund 2001), the metallicities range from 0.42 Z 
to 4.01 Z Q with a median of 1.69 Z Q . For the galaxies with 
measured Balmer decrements, 4 have metallicities closer to the 
L99 model with Z= 0.4Z Q and 47 have metallicities closer to 
the model with Z= 2Z Q . We test the differences between using 
the solar metallicity model and the model closer in metallicity. 
Assuming (B-R)q = 1.5, the maximum differences in Sr(i) and 
Wst are 0.06 and 15.9 Myr, respectively. Thus, we assume so- 
lar metallcity hereafter, which has no qualitative effects on our 
conclusions. 

Our data are insufficient to distinguish among different val- 
ues of (B-R)o. However, optical-to-near-infrared colors will 
aid in the process. For example, a galaxy like NGC 27 19 A, 
with measured EW(Ho;) = 100 A and (B-R) s = 0.56 has equally 
valid solutions for steep-IMF bursts of (among others) Sr(0 = 
0.54 and W = 23 Myr if (B-R) = 1.64 and s R (t) = 0.3, 
W = 12 Myr if (B-R) Q = 0.86. With a Salpeter-IMF burst 
on the same pre-existing blue galaxy, the solution differs even 
more significantly in age, with tburst = 329 Myr and Sr(i) = 0.48. 
Using Bruzual & Chariot (1996) Salpeter IMF t models for 
the pre-existing galaxies, with r = 1 Gyr and r = 100 Gyr (i.e., 
a constant star formation rate), respectively, we predict dif- 
ferences in the burst+pre-existing galaxy V — K colors of 0.17 
magnitudes between the two solutions with steep IMFs, and 
0.29 magnitudes between the steep IMF solution on the red 
pre-existing population and the Salpeter-IMF solution on the 
blue pre-existing population, with the redder V -H colors for 
the t = 1 Gyr Salpeter pre-existing population model and the 
steep-IMF burst. Thus, near-infrared colors help to distinguish 
among star formation histories and IMFs. 

Without near-infrared data, we must make assumptions about 
the properties of the pre-existing stellar populations. Fig. 10 
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shows the range of central-aperture colors in our sample. We 
segregate the sample into galaxies with substantial ongoing star 
formation (left, EW(Ha) > 10 A) and galaxies with little ongo- 
ing star formation (right, EW(Ha) < 10 A). We expect some 
contributions from triggered star formation even to the right 
side of the figure because timescales for blue B-R colors are 
substantially longer than timescales for Ha emission. Despite 
the contribution from triggered bursts, the outer edge of the dis- 
tribution (from interlopers and/or galaxies with no triggered star 
formation) provides an upper limit to the appropriate value of 
(B-R\). For illustrative purposes, we choose a fiducial constant 
value (B-R)o =1.5 (vertical line in Fig. 10b). 

Fig. 1 1 shows the 78 pair galaxies with measurable Balmer 
decrements (hence relatively strong ongoing star formation) on 
the continuous star formation model contours with (B-R)q = 
1.5 for three IMFs and three reddening corrections. The points 
fall in the expected range for the steep and cutoff Salpeter IMF 
slopes, based on all but the "traditional" reddening correction, 
which over-corrects continuum colors. For comparison, Fig. 12 
shows the contours and data with (B-R)o = 0.9. 

Although the absolute positions of the points on the contours 
are extremely uncertain, for a given (B-R)o the relative po- 
sitions of the galaxies in the (B-S) s /EW(Ha) meas (t) plane are 
more reliable. To the extent that a constant (B-R)o is valid, one 
conclusion from Fig. 1 1 is that the galaxies in our sample with 
larger EW(Ha) meas(0 have stronger and younger bursts of star 
formation. The more complex models that follow confirm this 
interpretation. 



3.4.1. Constraints on the IMF 

For all reddening corrections, most of the data fail to over- 
lap the Salpter model contours in Figs. 11 and 12. Most of 
the data that are bluer than the assumed (B-R)q fall below 
the constant-(Z?-/?)o Salpeter contours, corresponding to ages 
» 10 9 years. Timing arguments alone rule out this possibility 
for the majority of the data; thus, if the Salpeter slope is cor- 
rect, one or more of the other assumptions are invalid. Differ- 
ent constant (galaxy-independent) (B-R)q values do not affect 
the results; changes in (B-R)q only shift contours horizontally. 
Different star formation histories also fail to solve the problem 
convincingly. Replacing the continuous star formation mod- 
els with short (or instantaneous) bursts yields ages of <C 10 7 
years for all of the bursts. Although our data do not strictly 
rule out the instantaneous-burst scenario, interaction timescales 
are > 300 Myr, or ^50 times the duration of this burst; only a 
tiny fraction of our sample would then be bursting. This argu- 
ment does not rule out bursty star formation occurring quasi- 
continuously. In a time-averaged sense, a bursty star formation 
history mimics a steeper IMF — in the time between bursts, the 
highest-mass stars are preferentially removed. Thus, our data 
do not strictly distinguish between a steeper IMF and a shal- 
lower IMF with bursty star formation. 

For simpler (non-bursty) star formation histories, our mod- 
els and observations suggest either (1) a burst IMF with a 
high-mass slope steeper than the Salpeter IMF or, (2) bursts 
of star formation with a Salpeter- type IMF that produce smaller 
changes in the galaxy color. In other words, for Salpeter-IMF 
bursts to work, hence to lie on contours like those in Figs. 1 1 
and 12, the (B-R)o color must shift with each galaxy — its 
color must be determined not by the burst but primarily by the 
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pre-existing galaxy. 

Near-infrared colors will help to distinguish between these 
possibilities in the future. In addition, a different line of ar- 
gument, available to us now, results from the anti-correlation 
reported in BGK between EW(Ha) and pair separation on the 
sky, AD. In the next section, we adopt different distributions 
of (B-R)o from the literature and investigate the relationship 
between the burst IMFs and the possible origins of the AD - 
EW(Ha) correlation. 

4. THE ORIGINS OF THE AD - (B-R) s - EW(Ha) 
DISTRIBUTION 

The primary result of the spectroscopic study of BGK is the 
dependence of EW(Ha) and the other emission lines on the pair 
separation on the sky, AD. The equivalent widths also corre- 
late with velocity separation, AV. BGK find that galaxies in 
pairs with smaller spatial separations have a much larger range 
of EW(Hct), whereas galaxies with larger separations (up to 
~50 h" 1 kpc) almost always have equivalent widths < 30 A. 
BGK show that the correlation can be explained by a scenario 
in which a burst of star formation begins at a close pass, con- 
tinuing and aging as the galaxies move apart. In this picture, 
the dynamical timescales match star formation timescales (for 
the steep IMF). The picture is also broadly consistent with the 
hydrodynamical simulations of Mihos & Hernquist (1996); Tis- 
sera et al. (2002) for progenitor galaxies with shallow central 
potentials (no bulges). The biggest uncertainty in the BGK pic- 
ture is the contribution to EW(Hct) from reddening and from 
old stellar populations in the centers of the galaxies. Our anal- 
ysis of B and R photometry provides new constraints on this 
problem. 

Fig. 13 shows the steep IMF model (B-R)o =1.5 contours 



along with our data corrected according to the "Calzetti" red- 
dening law and segregated based on AD. We show only galax- 
ies with AD < 20 h" 1 kpc (filled circles) or AD > 30 h" 1 kpc 
(open circles). Although the results depend on the IMF and 
the reddening correction, the plot suggests that for the single- 
(B-R)q model, the bursts in the galaxies with small separations 
and large EW(Hct) are typically both younger and stronger than 
the bursts in the other galaxies, including the galaxies with 
larger separations on the sky. For a constant burst strength 
[sr(i) or Sioo], the galaxies in pairs with smaller AD tend to 
have younger bursts. 

In Sec. 3.4, we show that our model fails for Salpeter-IMF 
bursts if the progenitor galaxies have a very narrow (B - R)o 
range; the Salpeter model contours do not occupy the appropri- 
ate region in (B-R) s - EW(Hct) space. These arguments do not, 
however, rule out Salpeter bursts superposed on galaxies with 
a wide range of (B-R)q. Here, we construct a simple Monte 
Carlo simulation to test whether Salpeter- or steep-IMF bursts 
and correlations between AD and Sr(i) and/or AD and tburst 
can give rise to the appropriate observed ranges of EW(Ha) 
and (B-R) s and to the observed AD - EW(Ha) correlation 
[and relative lack of AD - (B-R) s correlation]. We restrict the 
comparison to the 72 galaxies with EW(Ha) > 15 A and mea- 
sured Balmer decrements, also applying this EW(Hct) limit to 
the simulated data. 

For the simulation, we draw 72 values of (B - R)q from the 
bulge colors of Peletier & Balcells (1996). We then draw a 
value of AD from a uniform distribution as indicated by the 
data; this uniform distribution represents an excess of very tight 
pairs over a random spatial distribution. For each tested max- 
imum burst age and strength, t max and s R . max , respectively, we 
explore the effects of correlations between pair separation, AD, 
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physical separation of the pair (or compact group). For a given 



FIG. 1 1 . — Pair data and model contours of constant sr(i) and t (see Fig. 3). 
Our model assumptions are described in Table 2, except we vary the IMF in 
the three columns: (left column) L99 steep IMF, (middle column) L99 Salpeter 
IMF, and (right column) L99 Cutoff Salpeter IMF. We correct the data for ex- 
tinction according to three prescriptions for dust (see Sec. 3.2.1): (top row) 
no correction (except for Galactic extinction), (middle row) "traditional" cor- 
rection, and (bottom row) "Calzetti" correction. We omit contour labels; they 
are the same as Fig. 3: (solid, left to right) sr(i) = 1,0.8,0.6,0.4,0.2,0.01 and 
(dotted, bottom to top) t = 10 9 , 10 85 , 10 8 , 10 7 5 , 10 7 , 10 65 , and 10 6 years. 
The lines on the right side of the plot show the median errors for the galaxies 
with measured Balmer decrements. For dust-corrected points, errors include 
propagated errors in the Balmer decrement. We plot the 112 points with no 
measured Balmer decrement on the top line only (open circles). 
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Fig. 12.— Same as Fig. 11 except (B-R)o = 0.9. 



and s R (t) and/or tburst- We test correlations that are "envelopes", 
appropriate because AD is only a constraint on the minimum 



AD, we draw s R (t) from to s R „ 
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AD 



This pro- 



4 (50 kpc) 

cedure results in a range of s R (t) values from to s R , max , with 
the maximum possible only at AD = for nonzero c s . Similarly, 



l-ct; 



AD 



tot rT 



we draw tburst from t max -(t max - 6 Myr)x - ^ (50kpc) 

resulting in a range of tburst values from 6 Myr to t max Myr, with 
the minimum possible only at AD = 0. The parameters c s and 
c t are the slopes of the "envelopes" and measure the strengths 
of the AD - s R (t) and AD - tburst correlations, respectively. We 
illustrate these "envelope" correlations in Fig. 14. We compute 
the "measured" EW(Hct) and (B-R) s directly from the selected 
(B - R)o, s R (t), and tburst using the L99 models with either the 
steep or Salpeter IMF. 

For several combinations of t max , s Rjmax , c s and, c t , we repeat 
this procedure 1000 times and examine the average. We ex- 
plore the ability of this procedure to reproduce values close to 
the averages, (EW(Ha)} and {(B-R) s ), of the data, corrected 
for extinction via the Calzetti, Kinney, & Storchi-Bergmann 
(1994) prescription (but not for nebular emission in this case); 
these averages are (EW(Hce)) = 61 A and ((B—R) s ) = 0.71, re- 
spectively. We also explore the relative strengths of the cor- 
relations we produced with this procedure using the Spearman 
rank probability as our measure of correlation. The correla- 
tion strengths for the data are somewhat ambiguous. With- 
out reddening correction or correction for nebular line emis- 
sion, the AD - EW(Ha) and (B-R) s correlations for these 
72 galaxies have strengths of P SR = 0.0077 and 0.033, respec- 
tively. With the "Calzetti" correction, which introduces scat- 
ter through imprecision but also renders (B-R)^ and EW(Ha) 
more accurate, the strengths of the Spearman rank correla- 
tions diminish to Ps R = 0.059 and 0.165, respectively. Be- 
cause the appropriate values of Ps R are somewhat unclear, we 
require only that the average Spearman rank probabilities sat- 
isfy 0.006 < (P SR ) < 0.06 for the simulated AD - EW(Ha) 
distribution, that 0.02 < (P SR ) < 0.2 for the simulated AD - 
(B - R) s distribution, and that the simulated AD - EW(Ha) 
relation is, on average, substantially more significant than the 
AD - (B-/?) s distribution averaged over the 1000 simulations. 

Fig. 15a shows (EW(Ha)) and ((B-R) s ) resulting from the 
simulations that satisfy the requirements on (Psr) for the AD - 
EW(Hct) and AD - (B-R) s distributions. We test only values 
with t max < 900 Myr. For all the tested Salpeter IMF bursts 
with the Peletier & Balcells (1996) distribution of (B-R) 
colors, the simulations fail to match the observed (EW(Ha)} 
and ((B-R) s ). The most pressing problem is that the initial 
bulge colors are too red; the EW(Hct) distribution requires weak 
Salpeter bursts for a match while the (B-R)^ constraint requires 
strong bursts. The combination that comes closest to the mea- 
sured range is that of strong (s R max = 0.99) but, on average, 
older bursts (t max = 900 Myr). No simple adjustments to the 
input AD - EW(Ha) or (B-R) s correlations will suffice for 
tmax < 900 Myr. If the Peletier & Balcells (1996) colors indi- 
cate the appropriate distribution of (B-R)o, the Salpeter bursts 
cannot reproduce the observations. 

Salpeter-IMF bursts on a bluer (B-R) Q distribution can meet 
with more success. We adopt the distribution of half-light col- 
ors from the Nearby Field Galaxy Survey of Jansen et al. (2000, 
NFGS), restricting to M B < -17 (with magnitudes as published 
but adjusted for Ho = 70 km s" 1 Mpc -1 ). Because colors are 
known to depend on M B , we ensure that the color distribu- 
tion drawn from the NFGS is appropriate for a sample with 
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FIG. 13. — Similar to the panel in Fig. 11 with the "Calzetti" dust correction 
and the steep IMF. The points are segregated based on pair separation, AD; 
rilled circles are galaxies in pairs with AD < 20 tr 1 kpc. Open circles have 
AD > 30 IT 1 kpc. We omit the galaxies without measured Balmer decrements. 
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FIG. 14. — An illustration of the simulated "envelope" correlations we draw 
from to construct the 72-galaxy simulated data sets that we compare with the 
real data. For this figure, we use c t = 4/9 and c s = 2/3 . The values we draw 
for the simulation are uniformly distributed in AD; hence they are not drawn 
uniformly from the shaded region, but uniformly from the x-axis. 



the M B distribution of our pairs. To match the distributions, we 
draw M B from the pair distribution; we then assign the galaxy 
a (B-R)o that equals the color of an NFGS galaxy selected 
randomly within the ±10 galaxies closest in M B . To take the 
most conservative approach, we correct the (B - R) Q colors us- 
ing Balmer decrements drawn from the pair sample, to match 
the E(B-R) in the pair bursts. This approach is conservative 
as it is appropriate only if the reddening of the pre-interaction 
galaxies equals the reddening of the post-burst galaxies. 

Fig. 15b shows the results from this algorithm. The Salpeter 
simulations with strong SR iiriax and large t max reach the range 
of interest for some simulations with t max = 700 or 900 Myr 
and SR max = 0.8. This maximum timescale is much larger 
that the expected dynamical timescales for ~50 kpc of sep- 
aration on the sky. They may reflect some contamination in 
the EW(Ha) > 15 A pair population from galaxies with older 
bursts or no triggered bursts, but the Salpeter IMF also prob- 
ably does not result in blue enough colors for the right range 
of EW(Hct). We note that the prescription of correcting the 
(B-R)o colors based on the measured burst Balmer decrements 
in our pair sample is probably an overestimate of the reddening 
for the pre-existing stellar population. In addition, the NFGS 
does not explicitly eliminate galaxies in pairs and therefore al- 
ready contains some triggered bursts. Finally, the (B-R)o color 
distribution we seek should include only the central few kpc of 
the galaxies, whereas the NFGS colors refer to the colors within 
the effective radius, frequently larger than the extraction area of 
our spectra; there is more blue disk contamination in the NFGS 
colors. Thus, only a conservatively blue (B-R)o color distribu- 
tion reproduces the observations with the Salpeter IMFs, with 
what are probably unrealistically long timescales. 

For the steep IMF, Fig. 15 shows that several combinations of 
parameters yield average EW(Hct) and (B-R) s values in the ap- 
propriate range, for either (B-R)o distribution. The best-fitting 



solutions we find all have s max > 0.5 and most have s max > 0.8. 
The solutions work well for a range of c s and c t values; all re- 
quire relatively strong AD - Sr(i) correlations, with c s > 0.33. 
However, the resulting timescales, with t max = 50 - 100 Myr, 
are unrealistically short. We therefore speculate that an IMF 
steeper than Salpeter, but shallower than the "steep" IMF we 
employ here, would probably provide the best match to the data 
with the most realistic timescales and assumptions. 

With the Peletier & Balcells (1996) bulge colors for the (B- 
R)o distribution, the best-fitting solutions all have AD - (B- 
R) s and AD - EW(Ha) correlations (with c s > 1 /3 and c t > 
1/3 for the solutions shown). However, for the bluer NFGS 
(B-R)q colors, 1 of the 13 indicated "steep IMF" solutions that 
fall within the indicated region has c t = 0. Thus, a AD - tburst 
correlation is not strictly required to reproduce the data. 

The Monte Carlo simulation is not a rigorous attempt to 
reproduce the data but a simple exploration of the different 
combinations of continuous bursts of star formation which can 
result in the appropriate AD - EW(Hct) - (B-R) s distribu- 
tion. We conclude that the result depends on the appropri- 
ate distribution of (B-R)o for the particular sample. If the 
Peletier & Balcells (1996) bulge colors are appropriate, the 
models require a high-mass IMF that is substantially steeper 
than Salpeter to explain the observed colors, equivalent widths, 
and AD - EW(Hct) correlation using simple star formation his- 
tories. Adopting much bluer progenitor colors does allow the 
Salpeter-IMF bursts to fit our model, but the required colors 
may be too blue to be realistic. Our data indicate a burst IMF 
steeper than Salpeter, but probably not as steep as a = -3.3. 

Regardless of IMF, the models explain the observed AD - 
EW(Hct) correlation only if there is a strong correlation be- 
tween pair separation on the sky and the fraction of /?-band flux 
in the spectroscopic aperture that originates from the new burst 
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<(B-R) s > <(B-R) S > 

FIG. 15. — Monte Carlo results using the (B-R) distribution of Peletier & Balcells (1996) (left), and from the NFGS (B-R)o distribution prescription (right) 
for which we draw (B—K)q colors based on Mb (see text). Solid points are Salpeter models; open circles show results using the steep IMF. The box indicates the 
region surrounding the sample average measurements (x), (EW(Ha)} = 61 A (±10 A) and ((B-R) s ) = 0.71 (±0.1), both corrected for reddening via the "Calzetti" 
reddening prescription. The points show Monte Carlo results for SR max = 0.2,0.5,0.8, and 0.99 (generally left to right), and t raax = 50, 100,200,300,500,700, and 
900 Myr (generally top to bottom), omitting combinations that result in very small (EW(Hct)). We add t max = 75 Myr for the Peletier & Balcells (1996) colors and 
steep IMF, to illustrate that the results coincide directly with the data. Considering values of c s and c t from to 1 (in increments of 1/9), we show only results that 
satisfy the Spearman rank correlation properties described in the text. 



(i.e., a strong AD - s R (t) correlation). Most, but not all, of the 
working models also require a strong AD - tb urs t correlation. 
The models which do not require the AD - tburst correlation 
represent a departure from the models of BGK, which explain 
the AD - EW(Ha) correlation with a AD - tburst correlation 
only. However, as we note, a substantial AD - tburst relation- 
ship is still consistent with the data (and predicted by most of 
the models we explore). If the correlation still plays some role, 
then the interpretation of BGK still holds qualitatively, but with 
somewhat shorter timescales than BGK discuss (e.g., by fac- 
tors of 1.5 - 3). We require additional data (e.g., near-infrared) 
to verify the relative roles of the AD - tburst and AD - s R (t) 
correlations. 

5. CONDITIONS THAT AFFECT THE STRENGTHS AND AGES 
OF THE BURSTS 

Although we lack direct constraints on (B-R)o we briefly 
investigate the relationships between third parameters and star 
formation more rigorously by inverting EW(Ha) meas (t) and (B- 
R) s to measure the burst age, tburst, and strength, s R (t). Several 
observational problems complicate this process. First, not every 
galaxy in the (_B-/?) s /EW(HaO me as(t) plane has a valid solution 
for every model. We omit galaxies with ages > 10 9 years. For 
the two galaxies with nominal s R (t) > 1 we define s R (t) = 0.999 
and estimate tburst by increasing the measured color until a valid 
solution exists; this procedure corresponds to moving the points 
to the right on Fig. 11, into the edge of the region of valid so- 
lutions. Second, uncertainties in (B-R)o plague this analysis; 
we focus primarily on models with constant (B-R)o and mod- 
els with a linear dependence of (B - R)q on the log of the rota- 
tion speed of the galaxy. The final observational problem re- 
sults from incompleteness in the sample of galaxies with mea- 



sured B aimer decrements (hence accurate corrected colors and 
equivalent widths); this sub-sample becomes incomplete for 
EW(Ha) meas(0 ^ 20 A (see Fig. 6). As a result of this incom- 
pleteness, the data are not sampled uniformly in tburst or s R (t). 
We explore possible correlations between third parameters and 
tburst and/or s R (t) after truncating the sample in accordance with 
this approximate incompleteness limit. There are many accept- 
able combinations of cutoff s R (t) and tburst that result in different 
sub-sample sizes. We explore two sub-samples: (1) s R (t) > 0.3, 
tburst < 10 8 years, and (2) s R (t) > 0.4, t bU rst < 10 82 years. Table 5 
lists the number of galaxies with valid solutions in the subsets 
we consider, including the full sample, the full sample restricted 
to galaxies with measured Balmer decrements, and these com- 
plete sub-samples. The table notes additional prescriptions for 
(B-R)o that we describe in Sec. 5.2, and indicates the number 
of galaxies that have measured velocity widths, 2 - 

5.1. s R (t), fburst, and the Orbit Parameters 

Fig. 16 shows the burst strength and age solutions assuming a 
constant (B-R)o = 1.5. We show all the galaxies with valid so- 
lutions with the steep IMF and the "Calzetti" reddening correc- 
tion. We plot burst s R (t) and age, tburst, as functions of AD and 
AV. The filled circles are the galaxies with measured Balmer 
decrements within the complete sub-sample restricted to galax- 
ies with s R (t) > 0.4 and tburst < 10 8 ' 2 years. The open circles 
indicate galaxies with solutions outside of the complete range 
and/or galaxies with no measured Balmer decrement. The typ- 
ical formal error in s R (t) is <~ 0.2 for 0.2 < s R (t) < 0.6 and is 
somewhat larger for other strengths. Errors in burst age are 
correlated with the burst strengths of the galaxies. For most 
galaxies, errors are ^0.2 - 0.3 dex; however, for galaxies with 
Sr(0 < 0.2 the age errors can span several dex (see Fig. 1 1). 
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Table 4 

Number of Galaxies with Valid Model Solutions for Subsets of the Data 



Prescriptions Number with Valid Solutions 









With 


Complete 2 


Complete 2 






C„11 

rull 


Measured 


(tburst < 10 yr) 


(tburst < 10 yr) 


(B-R) 


Reddening 


Sample 


T T /T T O 

Ha/Hfj 


(s R (t) > 0.4) 


(s R (t) > 0.3) 


(B-R) = 1.5 


"Calzetti" 


131 


77 


43 


51 


(B-R) =l.5+ 


"Calzetti" 


68 


35 


12 


16 


1.5[logV 2 c 2 -2.5] b 












CB -/?)„= 1.5+ 


"Calzetti" 


74 


39 


15 


20 


0.75[logV 2 c 2 -2.5] b 












(B-R) =1.5 


none 


114 


75 


17 


28 


(B-R) =1.5+ 


none 


57 


31 


1 


5 


1.5[logy 2 c 2 -2.5] b 












(B-/?)o=1.5+ 


none 


62 


37 


3 


8 


0.75[logV 2 c 2 -2.5] b 












Prescriptions 


Number with Valid Solutions and Measured V 2 2 


(B-R) = 1.5 


"Calzetti" 


72 


39 


17 


21 


(B-R){) =1.5 


none 


61 


38 


6 


11 



"Subsample in a "complete" region of tb urst - SR(t) space that includes only galaxies with measured Balmer decrements. 
*A11 galaxies in these categories have a measured V| 2 . 



The only visually evident correlation among the 43 solid 
points in Fig. 16 (complete sample) is the absence of very 
strong bursts with large separations. However, Spearman rank 
tests confirm both a tendency for galaxies with smaller separa- 
tions to have younger burst ages (Psr = 0.030) and for galaxies 
with smaller separations to have stronger bursts (Psr = 0.033). 
The other parameter combinations show no correlations, with 
P SR = 0.17 and 0.76 for the t - AV, and s R (t) - AV plots, re- 
spectively. If we instead adopt tburst = 10 8 years and s R (t) = 0.3 
as our completeness limits, the significance of the AD - tburst 
correlation increases to Psr = 0.0019; the significance of the 
AD - Sr(i) correlation also increases slightly, to Psr = 0.012. 
With the inclusion of the galaxies that lack measured Balmer 
decrements or that fall in the incomplete parts of the sample 
(open circles), a significant correlations persists between AD 
and s R (t) (Psr = 0.038) while the significance of the AD - tburst 
dependence drops to Psr = 0.12. 

5.2. (B-R)o, sr(i.,), /burst, and the Intrinsic Properties of the 
Progenitor Galaxies 

Sec. 4 examines general trends in the star formation timescales 
and orbital parameters of galaxies in pairs. Uncertainties in the 
properties of the progenitor galaxies, in the IMF of the triggered 
bursts, and in the reddening corrections complicate the inter- 
pretation of these apparent trends. Even if a theoretically tight 
relationship exists between the triggered star formation prop- 
erties and the positions of the galaxies in their orbits, we still 
expect variations in the progenitor galaxies (mass, orbit, envi- 
ronment, star formation history, IMF). These variations add to 
scatter introduced by interlopers (~10%; see BGK), by galax- 
ies that have not had a close pass yet, and by the accessibility 
of only two components of the true spatial AD and one of the 
true AV. 

One powerful tool for sorting out the causes of the large 
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FIG. 16. — Burst parameters (sR(t), and age, t, assuming (B-R)o = 1.5) vs. 
pair orbit parameters, AD and AV, for the galaxies with valid solutions only. 
Solid points represent the 43 galaxies in the complete sub-sample (excluding 
one point with AD = 76 IT 1 kpc from the left panels); they have measured 
Balmer decrements with burst ages t < 10 8 2 years and strengths SR(t) > 0.4. 
We plot the points without measured Balmer decrements and/or without valid 
solutions within the complete region as open circles. Within the complete sub- 
sample, there is a tendency for galaxies with smaller values of AD to have 
younger burst ages (Psr = 0.030) and stronger bursts (Psr = 0.033). 



scatter we find in Sec. 5.1 is the search for dependence of the 
amount and kind of interaction-triggered star formation on third 
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FIG . 17. — (B-R) s corrected for reddening and nebular lines (left) as a function of rotation curve velocity width and (right) R-band total magnitude for all galaxies 
in the SP sample. We segregate galaxies into groups with measured (solid points) or unmeasured (open points) Balmer decrements. The solid lines indicate the two 
prescriptions we adopt for (B-R)q as a function of 2 - 



parameters (e.g., Tissera et al. 2002). In this section, we explore 
the possible dependence of (B-R)q on intrinsic luminosity and 
rotation speed, and the relationships between burst parameters 
and the velocity widths of the galaxies. 

Galaxy colors depend strongly on luminosity, as exhibited by 
the different Tully-Fisher slopes observed for different band- 
passes (see e.g., Tully, Mould, & Aaronson 1982; Tully & 
Pierce 2000; Bell & de Jong 2001). However, estimates in 
the literature vary widely and depend on the particular veloc- 
ity width measure. Tully & Pierce (2000) find a difference of 
only ^0.3 magnitudes per dex; Bell & de Jong (2001) quote 
differences from 0.67 to 1.05 magnitudes per dex, depending 
on the reddening prescription. We concern ourselves only with 
the colors measured in the central (spectroscopic) apertures of 
the target galaxies. Thus, the extensive Tully-Fisher studies 
in the literature are of limited use. We rely on empirical es- 
timates of the probable range of variations in color with rota- 
tion speed. Fig. 17 shows (B-R) s , corrected for reddening with 
the "Calzetti" prescription, as a function of velocity width, for 
galaxies with well-measured rotation curves, and as a function 
of M R for all the SP galaxies. [See Barton et al. (2001) for 
the definition of the velocity width, which we measure from 
the major axis rotation curves.] The figure clearly shows the 
dependence of color on luminosity or rotation speed; the prob- 
able dependence of (B - R)o on luminosity and rotation speed 
follows. 

To test the effects of allowing (B-R)o to vary with veloc- 
ity width, we adopt the (B-R)q relations drawn in Fig. 17a. 
The lines correspond to (B-R)o = 1.5+1.5 [logV 22 -2.5] and 
(B-R) Q = 1.5 + 0.75 [log VS[ 2 - 2.5]. Using these prescriptions, 
we find that almost all of the correlations in Sec. 5.1 have been 



reduced in significance. The reduced sample size (from 43- 
51 to 12-20 galaxies; see Table 5) probably contributes to this 
reduction in significance. We note that a AD - tb urs t correla- 
tion persists, with P S r = 0.034, if we assume (B-R)q =1.5 + 
0.75 [logV 22 -2.5] and only accept solutions with tb ur .st < 10 8 
years and Sr(0 > 0.3. For the "complete" subsample with the 
same limits but with (B-R) Q = 1.5 + 1.5 [log V 22 - 2. 5], the cor- 
relation disappears, with Psr = 0.17. For the other "complete" 
subsample, these correlations also disappear, with Psr = 0.18 
and 0.33 for the two (B-R)q prescriptions. Curiously, we also 
note marginal correlations between AV and Sr(i) for various 
restricted sub-samples. 

The metallicities, emission line equivalent widths, and dust 
content of galaxies all depend on mass and/or luminosity (e.g., 
Tully, Mould, & Aaronson 1982; Skillman, Kennicutt, & Hodge 
1989; Heckman et al. 1998; Carter et al. 2001; Jansen, Franx, 
& Fabricant 2001). The required reddening correction thus de- 
pends on luminosity and rotation speed (e.g., Wang & Heck- 
man 1996; Tully et al. 1998; Sullivan et al. 2001; Hopkins et 
al. 2001). Fig. 18 illustrates this very significant dependence in 
the SP sample. The true fractional strengths of bursts of star 
formation in these galaxies may also depend on mass. Some 
trends with mass or luminosity appear to extend to the triggered 
star formation; in the BGK pair sample, the galaxies with the 
largest measured values of EW(Ha) are the lower-mass, lower- 
luminosity galaxies. Barton et al. (2001) find a shallow slope 
of the Tully-Fisher relation for these pairs. 

Here, we test for the dependence of the strengths and ages 
of triggered bursts of star formation on the masses (rotation 
speeds) and luminosities of the progenitor galaxies. Timescales 
for tidally-triggered star formation are much shorter than a 
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FIG. 18. — Balmer decrement as a function of R-band luminosity for (top) 
the 78 galaxies with measured Balmer decrements. The data show a tendency 
for more luminous galaxies to harbor more dust (P SR = 7.1 X 1CT 5 ). The his- 
togram shows the 112 galaxies with U/3 emission too weak to measure the 
Balmer decrement. These galaxies tend to be luminous. The histogram in 
Fig. 6 shows their EW(Ha) distribution. 



Hubble time. Thus, even if the frequency of interaction de- 
pends on mass, a large enough, unbiased set of galaxies in 
pairs should show no dependence of burst age on stable (non- 
transient) properties like progenitor mass. The Spearman rank 
test indicates, as expected, no dependence of burst age on rota- 
tion speed within the 17-galaxy complete sub-sample (see Ta- 
ble 5) when we assume a fixed (B-R)o =1.5 (although a cor- 
relation does exist between galaxy R-band absolute magnitude, 
M R , and tb urs t, with Psr = 0.011 or 0.031, depending on the 
complete sub-sample). We introduce very strong false corre- 
lations when we include galaxies in the incomplete region of 
Wst - Sr(0 space or galaxies without measured Balmer decre- 
ments; with the "Calzetti" extinction correction the significance 
is Psr = 0.0014. The complete sample with no dust corrections 
is too small to test for a significant correlation. The correla- 
tion including all the points in the no-dust-correction scenario 
is very strong (Psr = 3.1 x 10~ 6 ). Thus we conclude that restric- 
tion of the sample to the complete sub-sample and the correc- 
tion for dust are very important for eliminating artificial system- 
atic trends (see Fig. 18). Our prescription for changing (B-R)q 
as a function of V 2 2 also eliminates most of the artificial de- 
pendence of burst age on rotation velocity. 

Although correlations between galaxy mass and burst age are 
almost certainly unphysical, a correlation between mass and 
burst strength could be real. The low-mass galaxies clearly 
tend to have bluer colors and higher EW(Ha) values than the 
higher-mass galaxies. We investigate the solutions in detail, 
solving directly for the burst strengths and ages. We find some 
statistical evidence that the suggested correlation between ro- 
tation velocity and burst strength exists. Assuming a con- 
stant (B-R)q = 1.5, significant correlations exist both within 
the restricted, robust, complete sub-sample and for all of the 
data (including the galaxies with no measured Balmer decre- 



ment). With a "Calzetti" dust correction the Spearman rank 
probability is Psr = 0.035 for the complete sample (restricted to 
tburst < 10 8 2 years, s R (t) < 0.4; 17 galaxies), although it is only 
Psr = 0.086 for the other complete sub-sample. The complete 
no-dust-correction sample is too small to test for a significant 
correlation. However, including all the galaxies with valid solu- 
tions, measured V 2 2 > an d no dust correction, the Spearman rank 
probability is Psr = 4.0 x 10~ n . Because these are the prescrip- 
tions for which the unphysical burst age-luminosity and burst 
age-rotation speed correlations are strong, however, systematic 
problems introduced by incompleteness and the lack of knowl- 
edge of the true (B-R)o almost certainly affect interpretation of 
the results. 
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a function of rotation speed with (left) a "Calzetti" dust correction and (right) 
no dust correction. The small number of filled circles have measured Balmer 
decrements and valid age and strength solutions within a "complete region" 
(t < 10 8 ' 2 years, SR(t) > 0.4); open circles have no measured Balmer decre- 
ment and/or solutions outside the complete region. 



We test the robustness of the apparent burst strength-rotation 
speed correlation by assuming the prescriptions for (B-R)o 
with V 22 that we describe above. If we assume (B-R)o = 
1.5 + 0.75[logV 2 2 - 2.5] and "Calzetti" reddening, the com- 
plete sub-sample with tburst < 10 8,2 years and s R (t) < 0.4 yields 
a significant V 2 2 - sr(0 correlation (Psr = 0.0023). Restric- 
tion to the other sub-sample yields no correlation, however. 
For (B-R) = 1.5 + 1.5[logV5 2 -2.5] (see Fig. 19), the V° 22 - 
Sr(i) correlation remains significant for "Calzetti" dust (Psr = 
0.011); again, we observe no correlations within the second 
complete sub-sample. Inclusion of the points outside the com- 
plete sub-sample yields Psr = 0.49 (35 galaxies) and includ- 
ing the additional points without measured Balmer decrements 
yields Psr = 0.10 (for 68 galaxies). The correlation among all 
galaxies in Fig. 19b is significant, with Psr = 3.2 x 10~ 4 . We 
note that of all the combinations we consider in which (B - R)o 
depends on V 22 , artificial V 22 - tburst correlations are only 
found using the (B-R)q relation with a slope of 0.75 with no 
dust corrections. 

The relationship between V 2 2 and the true masses of galaxies 
is not straightforward, especially for interacting galaxies with 
centrally-concentrated star formation. Barton et al. (2001) find 
several examples of galaxies with truncated rotation curves that 
have anomalously small velocity widths. When analyzed from 
the perspective of another measure of intrinsic galaxy mass, 
Mr, the data do not present a particularly consistent picture. For 
example, with a constant (B-R)o = 1.5, we find no significant 
burst strength-luminosity correlations in the complete subsets 
of the data. When we restrict again to galaxies with measured 
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FIG. 20. — For the IRAS psc catalog detections, we plot a histogram of the 
log of the ratio of the total IRAS 60 ^tm flux within 3' of the galaxy to the R- 
band luminosity. We note that in some cases the IRAS flux includes flux from 
both galaxies. The shaded histogram includes only the galaxies with measured 
Balmer decrements (42 galaxies); the main histogram includes all galaxies 
with IRAS detections (120 galaxies). A K-S test indicates only marginally 
significant differences, at best, in the distributions, with P K s = 0.12. 



Vjj and assume either of the relations between (B-R) and 
Vj 2 we find a significant Sr(i) - Mr correlation for one of the 
complete sub-sample with tburst < 10 8 ' 2 and s R (t) > 0.4. In addi- 
tion, we note that there are unphysical Mr - tburst correlations, 
even among the complete sub-samples. We do note a strong 
tendency for galaxies with lower metallicities to have stronger 
bursts. Although we use solar metallicity for all of the bursts 
described here, our tests reveal that the use of lower metallicity 
models for the low metallicity galaxies and higher metallicity 
models for the higher metallicity galaxies would only enhance 
this correlation, as the better-suited models result in stronger 
bursts at low metallicity and weaker bursts at high metallicity. 
We plan to explore these relationships in more detail in the fu- 
ture. 

Our analysis, though far from conclusive, suggests that more 
luminous galaxies may experience weaker fractional bursts of 
star formation. However, as Heckman et al. (1998) point out 
and as Fig. 18 shows, more massive galaxies are more dust-rich. 
Strong dust-enshrouded bursts may hide among the luminous 
galaxies for which we are unable to measure a Balmer decre- 
ment. Hence the apparent conclusion (here and elsewhere) that 
less massive galaxies have higher luminosity-weighted star for- 
mation rates may be an artifact of inadequate correction for ex- 
tinction. 

The IRAS point source catalog provides a method of search- 
ing for strong starbursts embedded in the luminous galaxies. 
Because the IRAS beam is large, we sum all flux within 3' of 
the optical galaxy, including the flux for both galaxies in many 
close pairs. Thus, these ratios are actually upper limits. Fig. 20 
shows the distribution of the IRAS 60/im-flux-to-/?-band lu- 
minosity ratios for the IRAS-detected sample with measured 
Balmer decrements (37/78 detected; shaded histogram) and the 
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FIG. 21. — The relationship between Balmer decrement and the 60^im-to-K- 
band flux for IRAS-detected galaxies in pairs (see Fig. 20.) The bottom his- 
togram is the distribution of flux ratios for galaxies without measured Balmer 
decrements. 



sample as a whole (79/190 detected; open histogram). The dis- 
tributions are at best marginally different, with P K s = 0.12. We 
find no evidence for extreme hidden starbursts among the galax- 
ies with small EW(Ha) (hence without measured Balmer decre- 
ments). 
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FIG. 22. — Color profiles of the 6 galaxies with sr(i) > 0.8 (assuming 
(B—R)q = 1.5, the steep IMF, and the "Calzetti" reddening prescription). These 
galaxies and many others with strong central star formation show blue dips in 
or near their centers; few of the galaxies are accompanied by tidal tails, al- 
though some companions have tails. Fig. 23 identifies the galaxies and shows 
B-band images. 
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We use Fig. 21 to estimate the possible effects of hidden 
bursts of star formation on the sample of galaxies without mea- 
sured Balmer decrements. Based on the rest of the sample, 
their Balmer decrements are probably in the range of 2.88 - 
8, appropriate for the far-IR range in our sample. Thus, the 
"Calzetti" corrections to the B-R color and EW(Ha) may be as 
large as ^0.83 magnitudes and a factor of 2.5, respectively, for 
the galaxies with the greatest amounts of dust. Although they 
cannot affect the correlations among the galaxies in the sample 
with measured Balmer decrements, more accurate dust correc- 
tions could move the open points in Fig. 19 to much greater 
strengths, reducing the significance of the observed Vj 2 - sr(0 
correlation. Thus, although the IRAS data argue against ex- 
tremely strong hidden starbursts, they do not rule out subtle 
systematic dependences between galaxy luminosity and dust 
content that could give rise to a false correlation between burst 
strength and galaxy mass. 

In summary, we find some evidence that triggered star for- 
mation may be (fractionally) stronger in lower-mass galaxies. 
This result, if confirmed, would provide one possible mecha- 
nism for the stronger star formation that is typically observed 
in lower-mass galaxies. However, we find much less evidence 
for this trend when we replace V c 2 2 w i tn Mr as a measure of the 
intrinsic mass of a galaxy. In addition, known systematic trends 
in the dust content and the stellar populations as a function of 
rotation speed or luminosity of the progenitor galaxies may be 
affecting our results. More accurate measures of the dust cor- 
rections and pre-existing stellar population colors are required 
to explore these issues further. 

6. THE GALAXIES WITH STRONG CENTRAL STAR 
FORMATION 

As one product of the analyses we implement in this paper, 
we use the (B-fi) s /EW(Ha) meas (t) plane to identify populations 
of galaxies with particularly strong bursts or in particular stages 
of interaction. Fig. 22 shows the B-R color profiles of galax- 
ies with convincing evidence for strong tidally-triggered bursts. 
These six galaxies have s R (t) > 0.8 using the steep IMF and the 
"Calzetti" dust correction. We show images of the galaxies in 
Fig. 23. 

The profiles show blue dips in or near the centers; these 
dips are frequent in our sample. By-eye estimates indicate that 
22/185 (12%) of the luminous (M B < -17) non-AGN galaxies 
have blue central dips in the SP sample. In contrast, the NFGS 
contains only 8/147 (5%) of luminous galaxies with blue cen- 
tral dips. Because pairs were not explicitly excluded by the 
NFGS, some of the dips in the NFGS sample may also result 
from interactions. 

The SPH simulations of Mihos & Hernquist (1996) predict 
gas infall and centrally-concentrated star formation for bulge- 
less galaxies in the early stages of interactions (see the discus- 
sion in BGK). This centrally concentrated star formation has 
implications for the overall structural parameters of the galax- 
ies, especially at higher redshifts where low-surface-brightness 
features are difficult to detect. Barton et al. (2001) and Bar- 
ton & van Zee (2001) study the transient effects of central star 
formation on emission-line velocity widths and half light radii. 
They find a small number of outliers to the Tully-Fisher rela- 
tion in the BGK sample with very narrow linewidths for their 
luminosities. Follow-up HI observations with the VLA reveal 
that the radio linewidths of these galaxies are much broader; the 
narrow optical line widths result from poor sampling of only the 



central parts of the rotation curves. Barton & van Zee (2001) 
demonstrate that centrally concentrated star formation also ar- 
tificially reduces the half-light radii of galaxies. Thus measures 
of linewidths and radii may lead to underestimates of the intrin- 
sic masses and sizes of star-forming objects. 

Although they exhibit triggered star formation and unusual, 
lumpy morphologies, the target galaxies in Fig. 23 have no long 
tidal tails. At best, 2/6 have very stubby possible tails. Numer- 
ical simulations demonstrate that tidal tails are transient and 
occur only in prograde interactions (Toomre & Toomre 1972). 
Thus, studies that identify interactions on the basis of the pres- 
ence or absence of tidal tails are incomplete. 

Our analyses also provide a method of identifying galax- 
ies that have already undergone a close pass and are slow- 
ing down for a second pass or merger. We select galaxies 
that have been experiencing a long-lived but ongoing burst of 
star formation and are near apogalacticon by choosing galax- 
ies with blue colors [here, (B-R) s c < 0.75], moderate EW(Ha) 
[20 A< EW(Ha) < 40 A] indicative of a substantial build-up of 
R-band continuum, and small velocity separations (AV < 150 
km s" 1 ). The morphologies of the chosen pairs suggest that this 
method is nearly 100% successful at identifying galaxies that 
have already experienced a close pass. We show B-band im- 
ages of the pairs in Fig. 24; one of the galaxies overlaps with 
the strong -burst sample of Fig. 23. Of the 7 galaxies we se- 
lect with these criteria, which are members of 7 different pairs, 
each pair shows one or more relatively unambiguous indication 
of a previous close pass (e.g., lumpy, unusual morphology or a 
tidal tail either in the blue galaxy itself or in the neighbor). A 
large sample of galaxy pairs near apogalacticon chosen in this 
manner would provide kinematic constraints on the dark matter 
content of galaxies (e.g., Charlton & Salpeter 1991). 

7. CONCLUSION 

We combine B and R photometry with longslit and cen- 
tral spectroscopy to explore the strengths and ages of tidally- 
triggered bursts of star formation in galaxies in pairs. We con- 
struct a two-population description of the pre-existing and trig- 
gered stellar populations near the centers of the galaxies. Our 
primary conclusions are: 

1 . The most realistic starburst models that describe the 
data statistically include continuous star formation, a 
steep (or truncated Salpeter) IMF, and the reddening 
correction of Calzetti, Kinney, & Storchi-Bergmann 
(1994). Consistency with the Salpeter IMF requires a 
very blue distribution of progenitor-galaxy colors. 

2. In our picture, the models require strong triggered 
bursts in some of the pair galaxies; these bursts 
constitute >50% of the central R-band light. This star 
formation leads to a high incidence of galaxies with 
blue centers; tidal tails do not always accompany these 
central bursts of star formation. 

3. Our results suggest that the strengths, and most 
likely the ages, of triggered bursts of star formation 
depend on the galaxy separation on the sky. Thus, 
although they do not conclusively verify the model, 
our data are consistent with the picture of a burst 
of star formation triggered by a close galaxy-galaxy 
pass that continues and ages as the galaxies move 
apart. In this picture, the strongest bursts of star 
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formation occur only in the tightest orbits, giving rise 
to the strength-separation correlations. The data also 
suggest a possible correlation between the strength 
of a triggered burst and the rotation speed (mass) 
of the progenitor galaxy, in the sense that low-mass 
galaxies experience stronger bursts of triggered star 
formation. This result supports the hypothesis that the 
evolution of galaxies is mass-dependent. Verifying 
these conclusions requires additional constraints on 
the colors and dust content of the pre-existing stellar 
populations. 



The colors of the pre-existing stellar populations and the 
reddening corrections are the dominant uncertainties affecting 
conclusion (3). More accurate characterization of the correla- 
tions between orbit parameters and star formation properties in 
tidally interacting galaxies requires independent constraints on 
the colors of the old stellar populations. Near-IR observations 
can provide more direct measurements of properties of the older 
stellar population; we plan to report soon on a comprehensive 
UB VRJHK imaging study of a subset of these galaxies in pairs. 



We thank Jan Kleyna, Norman Grogin, and Daniel Koranyi 
for useful software and help with the data reduction procedures. 
We thank Lisa Kewley for the use of her abundance modeling 
code, and for many helpful comments; we also thank Eric Bell, 
Liese van Zee, Rob Kennicutt, and Rolf Jansen for helpful com- 
ments relating to this work. We are grateful to Perry Berlind and 
Mike Calkins for assistance with the spectroscopic observations 
and to Susan Tokarz for assistance with data reduction. Support 
for E. B. G. was provided by NASA through Hubble Fellowship 
grant #HST-HF-01 135.01 awarded by the Space Telescope Sci- 
ence Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., for NASA, under con- 
tract NAS 5-26555. The Smithsonian Institution supports the 
research of M. J. G. and S. J. K. 



REFERENCES 

Abraham, R. G., Tanvir, N. R., Santiago, B. X., Ellis, R. S., Glazebrook, K., & 

van den Bergh, S. 1996, MNRAS, 279, 47 
Allende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 556, L63 
Andredakis, Y. C, & Sanders, R. H. 1994, MNRAS, 267, 283 
Barton, E. J., & van Zee, L. 2001, ApJ, 550, L35 

Barton, E. J., Geller, M. J., Bromley, B. C, van Zee, L., & Kenyon, S. J. 2001, 
AJ, 121, 625 

Barton, E. J., Geller, M. J., & Kenyon, S. J. 2000a, ApJ, 530, 660 
Barton, E. J., Kannappan, S. J., Kurtz, M. J., & Geller, M. J. 2000b, PASP, 
112,367 

Bell, E. R, & de Jong, R. S. 2000, MNRAS, 312, 497 
Bell, E. R, & de Jong, R. S. 2001, ApJ, 550, 212 
Bessell, M. S. 1990, PASP, 102, 1181 
Buser, G., & Kurucz, R. L. 1978, A&A, 70, 555 

Bruzual, A. G., & Chariot, S. in AAS CD-ROM Series, Vol. 7, Astrophysics 

on Disc (Washington: AAS) (GISSEL96) 
Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582 
Carollo, M. C. 1999, ApJ, 523, 566 

Carter, B. J., Fabricant, D. G., Geller, M. J., Kurtz, M. J., & McLean, B. 2001, 

ApJ, 559, 606 
Chariot, S., & Fall, S. M. 2000, ApJ, 539, 718 
Charlton, J. C, & Salpeter, E. E. 1991, ApJ, 375, 517 
Courteau, S. C. 1997, AJ, 114, 2402 

Diaferio, A., Kauffmann, G., Colberg, J. M., White, S. D. M. 1999, MNRAS, 
307, 537 

Donzelli, C. J., & Pastoriza, M. G. 1997, ApJS, 111,181 

Falco, E. E., Kurtz, M. J., Geller, M. J., Huchra, J. P., Peters, J., Berlind, P., 

Tokarz, S., & Elwell, B. 1999, PASP, 111, 438 
Heckman, T. M., Robert, C, Leitherer, C, Garnett, D. R., & van der Rydt, F. 

1998, ApJ, 503, 646 



Hopkins, A. M., Connolly, A. J., Haarsma, D. B., & Cram, L. E. 2001, ApJ, 
122, 288 

Jansen, R. A., Franx, M., Fabricant, D., & Caldwell, M. 2000, ApJS, 126, 271 
Jansen, R. A., Franx, M., & Fabricant, D. 2001, ApJ, 551, 825 
Jones, B., & Stein, W. A. 1989, AJ, 98. 1557 

Kauffmann, G, Colberg, J. M., Diaferio, A., & White, S. D. M. 1999, 

MNRAS, 303, 188 
Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 1999, 

MNRAS, 307, 529 
Keel, W. C. 1993, AJ, 106, 1771 
Keel, W. C. 1996, ApJS, 106, 27 
Kennicutt, R. C, & Keel, W. C. 1984, ApJ, 279, L5 
Kennicutt, R. C, Jr., Keel, W. C, van der Hulst, J. M., Hummel, E., & 

Roettiger, K. A. 1987, AJ, 95, 5 
Kennicutt, R. C, Jr., Tamblyn, P., & Congdon, C. W. 1994, ApJ, 435, 22 
Kewley, L. J., & Dopita, M. J. 2002, ApJS, in press (astro-ph/0206495) 
Koo, D. C, Bershady, M. A., Wirth, G. D., Stanford, S. A., & Majewski, S. R., 

1994, ApJ, 427, L9 

Landolt, A. U. 1992, AJ, 104, 340 

Larson, R. B., & Tinsley, B. M. 1978, ApJ, 219, 46 

Leitherer, C, et al. 1999, ApJS, 123, 3 

Liu, C. T., & Kennicutt, R. C, Jr. 1995a, ApJ, 450, 547 

Liu, C. T., & Kennicutt, R. C, Jr. 1995b, ApJS, 100, 325 

Madore, B. F. 1986, in Spectral Evolution of Galaxies, eds. C. Chiosi & A. 

Renzini (Reidel, Dordrecht), p. 97 
Mihos, J. C, & Hernquist, L. 1996, ApJ, 464, 641 
Miller, J. S., & Mathews, W. G. 1972, ApJ, 172, 593 
Miller, G. E., & Scalo, J. M. 1979, ApJS, 41, 513 
Peletier, R. F, & Balcells, M. 1996, AJ, 111, 2238 
Pfenniger, D., & Norman, C. 1990, ApJ, 363, 391 
Salpeter, E. E. 1955, ApJ, 121, 161 
Sekiguchi, K., & Wolstencroft, R. D. 1992, 255, 581 
Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 
Skillman, E. D., Kennicutt, R. C, & Hodge, P. W. 1989, ApJ, 347, 875 
Sullivan, M., et al. 2001, ApJ, 558, 72 

Thompson, R. I., Weymann, R. J., & Storrie-Lombardi, L. J. 2001, ApJ, 546, 
694 

Tissera, P. B., Dommguez-Tenreiro, R., Scannapieco, & Saiz, A. 2002, 

MNRAS, 333, 327 
Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 
Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 
Tully, R. B., Mould, J. R., & Aaronson, M. 1982, ApJ, 257, 527 
Tully, R. B., Pierce, M. J., Huang, J.-S., Saunders, W., Verheijen, M. A. W., & 

Witchalls, P. L. 1998, AJ, 115, 2264 
Tully, R. B., & Pierce, M. J. 2000, ApJ, 533, 744 

Veilleux, S., Kim, D.-C, Sanders, D. B., Mazzarella, J. M., & Soifer, B. T. 

1995, ApJS, 98, 171 

Wang, B., & Heckman, T. M. 1996, ApJ, 457, 645 

Zombeck, M. V. 1990, Handbook of Space Astronomy & Astrophysics, 

Second Edition (Cambridge: Cambridge University Press) 
Zwicky, F. & Kowal, C. T. (1968). Catalogue of Galaxies and Clusters of 

Galaxies (California Institute of Technology, Pasadena). Vol. 6 




FIG. 23. — B-band images of the galaxies in Fig. 22. We center the images on the galaxies with nominal Si$(t) > 1. Galaxies "A" - "E" are UGC 312, IC 2338, 
NGC 3991 -North, CGCG 132-062, and UGC 1201 1; they are labeled in Table 1. "D" refers to only the Northern galaxy in the close galaxy pair (middle right). 
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FIG. 24. — Pairs with one member selected based on SR(t), EW(Ha), and AV to be near apogalacticon. Although some faint features are difficult to see in this 
figure, most of the pairs show signs of tidal interaction and/or lumpy star formation. The selected galaxies, "G" - "M", are CGCG 059-050, CGCG 180-020, NGC 
3991-North, UGC 07085-West, CGCG 217-024, CGCG 132-063, and UGC 12265-South; they are labeled in Table 1. (Note that "D" in Fig. 23 is the same galaxy 
as "I" here and "L" is the companion of "E".) 



